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ABSTRACT 

The effect of finite ion Larmor radius on the develop- 
ment of Rayleigh-Taylor instability of a plasma bounded 
by a vacuum is studied. The macroscopic equations of mo- 
tions are used where the finite ion Larmor radius effect is 
incorporated through off-diagonal terms in the pressure 
tensor in the momentum fluid equations. The effect of the 
finite Larmor radius stabilization of the interchange mode 
is demonstrated. For disturbances propagating along the 
magnetic field it is found that the inclusion of gyro-viscosity 
has the effect of reducing the instability of plasma vacuum 
interface. 
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I. INTRODUCTION 

The importance of the effect of the finite Larmor radius (FLR) of the ions 
on plasma instabilities has been of considerable interest in the recent years. 

The effect of FLR on the gravitational instability of a plasma was considered by 
Rosenbluth et al 1 using Vlasov equation and the same results were recovered 
later by Roberts and Taylor 2 using macroscopic equations for a two dimensional 
plasma. Recently Kennel and Greene 3 have investigated the FLR effects using 
the appropriate scalings of the Vlasov equation and derived a set of asymptotic 
equations. Among several interesting results, they obtain the FLR stabilization 
of the interchange mode. 

This paper considers the stability of an incompressible plasma bounded by 
a vacuum and supported against gravity by a static magnetic field. We use here 
the macroscopic moment equations where the effect of finite ion Larmor radius 
is included through off-diagonal terms in the pressure tensor in the momentum 
fluid equations. The boundary conditions relevant to the problem at hand are 
derived and the dispersion relation is obtained. The interchange mode stabili- 
zation is demonstrated. When the disturbance is propagating along the magnetic 
field, the dispersion relation is analyzed on the assumption that gyroviscosity is 
small, i.e. NDkV/Onw, « |, where N is the particle density, ©the plasma temper- 
ature in energy units, k and n are the wave number and frequency of the dis- 
turbance respectively, p the density and u> c the ion-gyrofrequency. It is found 
that the inclusion of gyro-viscosity has the effect of reducing the instability of 
plasma vacuum boundary. 
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II. FORMULATION OF THE PROBLEM. 


We now consider a situation where an infinitely conducting plasma occupies 
the half space 0 < y <oo and is supported against gravity by a uniform magnetic 
field which we shall take along the z-axis of a system of cartesian co-ordinates 
(x, y, z). We shall assume that the medium is incompressible and infinitely 
conducting. The equations basic to our discussion are the equations of motion 

P — — - — — V 77 + — — (B ■ V) B + p g — V' P M • (1) 

dt 4 77 

and the generalized Ohm's law 

E + ivxB - _2 (V x B) x B + -L Vp = 0 (2) 

c 4 77Ne Ne e 

together with 

V • v = 0 and V * B = 0, (3) 

where v = p + B 2 /87t, g = (0, -g, 0) is the constant gravitational field, N is 
the electron number density, p e is the electron pressure, p is the total fluid 
pressure, p is the constant density, v is the fluid velocity, E the electric 
field, B the magnetic field and c the velocity of light. The electric and the 
magnetic fields satisfy Maxwell's equations in which we shall neglect the dis- 
placement current. The anisotropic part of the stress tensor, P M , arises because 
of the finite Larmor radius effects. If the magnetic field is along the z-axis, its 
various components are given by 4 


4 


P = -ov A , P - p vk , P ~ 0, 

XX r xy yy r xy zz 


( 4 ) 


P = P = ~2pvk , 

zx xz r yz 

P = P = 2pvk , 

zy yz ' zx 

P = P = 1. pv (A - A ), 

xy yz y H v xx yy y 


(5) 

(6) 
(7) 


where 


A = I [Vv + (Vv) tr ] - I V *v I, (8) 

" 2 3 

and I is the unit dyadic. The coefficient v occurring in the above relations is 
given by v =N©/ y0 6J c , where 0 = kT and = eB/Mc is the ion gyrofrequency. It 
is clear from the foregoing equations that these admit the static stationary 
solution 


V n - Vp = p g . 


(9) 


The condition of pressure balance in equilibrium at the interface then becomes 



8 77 



( 10 ) 


where B p and B v denote the strength of the magnetic fields in the plasma and 
in the vacuum respectively. 
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For small departures from the state of equilibrium, let the various quantities 
be denoted by v, B p + b, p + Sp. The equations governing these are readily ob- 
tained from Eqs. (1) - (4) together with Maxwell's equation BB/Bt = -cV x E. 

We thus obtain 


and 


Bv - 1 

— — = - V 7 7, + 

Bt 1 477 p 


(B„ 


= Vx (vxB)-c, 
Bt v p ; i 


V: 


• V)6-lv -p u , 

[(Vx b) x B p ] , 


(ID 

( 12 ) 


where 77 x = Sp / p + B p • b/ 4-np and c x = c/477Ne = B p /4 rrpco^ . 

We now assume the (t,z) dependence of the various quantities to be of the 
form f(x, y , z, t) = f ( x, y ) exp (n t + ik z z). Equations (11) and (12) then reduce 
to 


and 


_ - i ik B _ 

nv = — V 77, - _ V • P„ + z p . b 

1 P M 477 p 


nb = ik B v + ik B c, V x b . 

Z p Z p 1 


(13) 


(14) 


The gradient term in Eq. (13) can be eliminated by taking the curl of this equa- 
tion. Assuming that the temperature variations in the initial state are negligible, 
we find after some straight forward reductions that 


nV xv = ik — (V x V x v - 3k 2 v) + — V x b ■ 

2 z 4rrp 


ikB - 

p 


(15) 
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We can then eliminate the velocity field from Eq. (15) by using Eq. (14) and ob- 
tain the equation for perturbation in the magnetic field b : 

a k -3 curl 3 b + a k' 2 curl 2 b + a,k'* curl b + a n b = 0, (16) 

o z 1 z 1 z U 

where we have put 


and 



nco v 

O 2 



B 2 

co = vV. 2 , a 2 = k 2 A 2 , A 2 = E 

* 4tt p 


III. THE METHOD OF SOLUTION. 


The solution of Equation (16) can be written as 


E =L B1 - 

i=i 


where b* is governed by the equation 


V x 9 = a. b l 


(17) 

(18) 

(19) 

( 20 ) 

(21) 


( 22 ) 


(23) 
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and a. ’s are solutions of the cubic equation 


a 3 k ; 3 + a 2 K 2 a i + a i k z 1 a i + a 0 n °- 


To solve Equation (23) we further assume that the x-dependence of the per- 
turbation is of the form ik x x. We then obtain for the solution valid in the upper 
half plane (y >0): 


- K . y 

b 1 = C.p. e 1 , 

X 1 *1 * 


y 

b 1 = iC.q.e 1 , 


y i ‘i 


where 


=C.e 


\ 2 - k 2 - a 2 , k 2 = k 2 + k 2 


k k + a. k. k a. + k 

XZ T 11 x i 1 Z 1 

Pi = q 5 = 


k 2 - a 2 

Z 1 


k 2 - a 2 

Z 1 


and C are the three integration constants to be determined by the boundary 
conditions. Having obtained b the perturbation in the velocity field can be ob- 
tained from Equation (14) to be 


5 = L = rr £ I" 

z p 1 = 1 


i Q 2 a. 
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The perturbation in the total pressure tt x can be obtained by taking the z-component 
of Eq. (13) and making use of the foregoing results. We thus find 


V. - / 77 . = 

L— j 1 B 

1=1 Z P i=l 


1 3 

tL 


iQ 2 a. 

n 2 + ( n + - ) (n - i vk z a. ) 

k co. 

Z 1 


(30) 


Finally the perturbation in the vacuum magnetic field, b ( °) is the solution 
of equations V x b (0) =0 and V • b (0) = 0; it is given by 


g(0) _ [ e ^ ik x + e y k + e z i k z ] C 4 e ky (y < 0), (31) 

where C 4 is another integration constant to be determined by the boundary 
condition. 

IV. THE BOUNDARY CONDITIONS 

The dispersion relation can now be obtained by using the boundary conditions 
appropriate to the problem. These are: (1) At the interface y = 0, the normal 
component of the velocity must be compatible with the assumed form of the 
deformed boundary, i.e., 


i ( k x + k z + n t ) 

v = n e e 


(32) 


where e is a constant. 

(2) The normal component of the magnetic field must be continuous. The 
linearized form of this condition is 
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+ B -8N] = 0. 


(33) 


A [b • N 0 

where A [X] denotes the jump in the quantity X at the interface. Here N 0 is the 
unit normal to the undisplaced interface and 8N denotes its displacement. For 
the problem at hand 

N 0 = e y and SN = - e (e x ik x + e z ik z ). (34) 

The remaining conditions are obtained by integrating the equation of motion 
(1) across the interface. We thus obtain 


N A (tt) = N • A (B B - P) ■ 


(35) 


(3) If we multiply Eq. (35) scalarly with N we obtain the boundary condition for 
the continuity of the normal component of the stress: 


A(tt) = N • A (N • P). 


(36) 


The linearized form of this equation leads to the condition 

^ - £p% + P yy =Tt\ - JL b(°> , (y = 0), (37) 

where the superscripts ‘p’ an d *v’ on 7Tj refer to the plasma and the vacuum 
quantities and in writing Eq. (37) we have made use of the relation 37 tp /By = -pg 
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The other two conditions are obtained by taking the vector product of Eq. 
(35) with N . This leads to the conditions for the continuity of the tangential 
stress: 


(N • B) N x A(B) = N x A (N . P) (38) 

The linearized form of Eq. (38) leads to the conditions 

P xy = 0 (at y = 0) (39) 

and 

A [B p (b y + B p SN z )] =P zy (y = 0). (40) 

On making use of the foregoing boundary conditions we obtain the following 
characteristic value equations: 


L 


(n + /3 a. ) q. C. = nek B , 

v 'i'll Z p ’ 


(41) 


22 * Ci = ~ ikC 4 +k z e ( B p- B v)’ 


22 (° + P a i) ( k x p i + \ q i) - o> 


(42) 


(43) 


22 V( ~ n + /3a i ) (kx + k z Pi) C i = -y- (B P " Bv)(k C 4 - 1 k Z 6 B v> ’ 


(44) 


11 



( 45 ) 



where we have put /3 = ifi 2 /k z a> c . Equations (41)-(45) are five equations involv- 
ing the five unknown constants C it (i = 1, 2, 3, 4) and e. The dispersion relation 
is obtained by setting the determinant of their coefficients equal to zero. 


V. THE INTERCHANGE MODE. 

In view of the complexity of the general case we first consider the simple 
case of the interchange mode. We observe that in this case k z = 0, k x / 0 and 
Equations (13) and (14) reduce to 


and 



(46) 

(47) 


Equation (47) gives b = 0, since n / 0, while Eq. (46) leads to 



where we put v = (e x u + e y v + e z w). The z-component of Eq. (48) requires 
that w = 0. The other two equations together with the divergence condition yield 
the solutions 
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• „ -k x y + ik x x + nt 
u = - i C e x 


( 49 ) 


v = C e 


c x y + ik x X + n 


(50) 


and 


nC - k x y + ik x+nt 
- — 6 


(51) 


The constant C is determined by the requirement that at y = 0, v = ne and we get 


C = ne. 


(52) 


The boundary condition (37) now reduces to 


1 I ■ , 3u 

1 2 I x 3y 


eg at y = 0 . 


(53) 


On using the foregoing results in this equation we obtain the dispersion relation 


n 2 = g k - invk 2 . 

X X 


(54) 


This leads to the result 


n = ± (gk x ) 1/2 



1/2 


i vV? 

X 

2 


(55) 


We may observe here that when v - 0, n =± (gk x ) 2 /2 which agrees with the 
classical result that the presence of a magnetic field has no effect on the inter- 
change mode. However, when ^ i 0, we find from relation (55) that (a) the 


13 



instability sets in as overstability and (b) the growth rate of instability is re- 


duced due to the gyro-viscosity; and this effect is large for short wavelengths. 
In fact we can see from Eq. (55) that all wave numbers with k greater than a 
certain critical k t are stabilized by the gyro-viscosity, where k„ is given by 



and all wave numbers k < k„ are unstable. 


(56) 


VI. THE CASE k x = 0 AND k z / 0. 

We shall now drop the subscript zonk as it is no longer necessary. Equa- 
tions (41)-(45) now lead to the dispersion relation: 


(n + /3 dj) (n+/3a 2 ) (n+/3a 3 ) 


Qi 


(n+/3a x ) qi (n+/3a 2 )q 2 (n+/3a 3 )q 3 


n 


= 0, 


(57) 


where 


1 - V 


- rjCl 1 2 kg 


Qi =/3a i q. - 


mpv 


(1 - v) 


(n + P a.. ) p. , 


(58) 


Pi = + (n + /3a.) [ n - a. 


(59) 
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and 17 = B v /B > 1. In order to solve the dispersion relation (57), we now need 
the expressions for a. ’s, the roots of Eq. (24). As it seems quite difficult to 
find the general solution of Eq. (24), we shall restrict ourselves to the case 
when the gyro-viscosity effects are small and we shall assume that w, /n « 1. 
The roots of Eq. (24) are then determined by an iteration procedure and we find 
that correct to the first order terms in a> v /n, the roots are given by (a> v = vV 2 ): 


a 


1 


ika; 

c 


n 



ik cx> 

V 

2n 



3 

+ 

(1 + n 2 /Q 2 ) 


and 


-j inkw 
3 v 

2 n 2 +n 2 


a 3 = - 


2ink 


f 1 CO CO 

- 1 V c 

1 + - 


CO 


(60) 


(61) 


(62) 


The appearance of the term (l/to v ) in Eq. (62) leads to no difficulties as we de- 
mand that in the limit v - 0, the constant C 3 - 0. With a’s given by Eqs. (60) - 
(62) we can readily see that to the lowest significant order, we have 


p >'£ rv ’(‘-r sx )- 


(63) 





So; 2 

% 

= r 

1 

c 

+ 



- 

2 n 2 



r 

S to 2 


k 

1 

c 

1 

"r 


2 n 2 


= k, 

^2 

= 1 , p. 


(i + 2r 2 ) 


3 . 8na; c 
2 n 2 a 2 


(64) 


(65) 


( 66 ) 
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( 67 ) 



and e . is the unit tensor of third rank which is completely antisymmetric in 
all the indices. On using the relations (60) - (70) in equations (71) - (73) and 
carrying out all the expansions consistently to the lowest significant order in S , 
we obtain after some straightforward reductions for the dispersion relation: 
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(74) 


D 0 + — 
0 2 


n 


2 ( co 2 o co 2 

lD 0 a 2 +D 0 ( 2 r 2 _JL - + _JL a 

2 \ n 2 a 4 n 2 


+ _ (N - n 2 ) (N + 77 n 2 ) 


where 


co u oj . n 

+ — — (N + V n 2 ) N - — _ _ 

n a 2 n n r ._•> 


n r a 2 


= 0, 


D 0 = kg-n 2 - Q 2 (1 + t? 2 ) , 


N = kg + rjQ. 2 (1 - r, 2 ), 


(75) 


2nQ. 2 (1 -r]) 


and b = 1 - a . 


First of all we observe that in the limiting case when S - 0, i.e., the ab- 
sence of gyro-viscosity, Eq. (74) reduces to the classical result 


D 0 = kg - n 2 - fi 2 (1 + i 7 2 ) = 0 , 


(76) 


where we put a subscript on n to denote its value in the limit b - 0. In order to 
obtain the first order correction to the dispersion relation, we substitute n 0 
given by Eq. (76) into all the terms containing S in Eq. (74) and obtain the result 


n 2 = kg - fi 2 (1 + 7 ? 2 ) - y- n 0 f , 


= n 


2 

0 



(77) 
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where 


f _ i _ (^ + 1) 

2 tt (j) - 1 ) 

Note that f is positive only when 77 is larger than ?] c = (2 tt+ 1)/(2t7- 1) 1 .4 • 

If we consider a situation where v > V c and which is unstable in the absence 
of co v i.e., n Q > 0, then for this unstable mode in the presence of co v , n 2 > n 2 > 0 
i.e., the instability is reduced. If on the other hand we consider a situation 
which is stable in the absence of co v i.e., n 2 < 0, then the inclusion of gyro- 

effects in the equations of motion leads to damping of the oscillations. 

VII. CONCLUSIONS 

We find that for the interchange mode with the inclusion of gyro-viscosity 
(a) the instability sets in as overstability and (b) the growth rate of instability is 
reduced. For disturbances propagating along the magnetic field, for a high /3- 
plasma (77 > v c ) the instability is reduced due to the inclusion of gyro-viscosity. 
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